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OBSERVACIONES AL ALUMNO 

Al comienzo de cada texto encontrará lo siguiente. 

1. UNA PRESENTACION DE NUESTROS OBJETIVOS en la unidad, la cual puede usar para comprobar si ha completado el 
curso satisfactoriamente. 

2. UN DIAGRAMA ESTRUCTURAL que le muestra cómo están relacionadas las diferentes secciones de la unidad. El diagrama 
consta de: 
a) cuadros rojos que muestran la secuencia principal. 
b) cuadros negros punteados (a la izquierda del rojo) que le muestran lo que esperamos ya sepa antes de fijarse un estudio 

detallado. 

e) cuadros negros (a la derecha del rojo) que le ayudan a escoger temas que puede posponer u omitir si tiene poco tiempo. 

3. UN GLOSARIO de las definiciones, notaciones, etc., usadas en la unidad como referencia inmediata. 

4. UNA BIBLIOGRAFIA de los libros que puede consultar cuando haya acabado el trabajo de la unidad. 

5. CONCLUSIONES Y RESUMEN del trabajo realizado al final de cada texto. 


NUMERACION E INDICADORES 

Dentro de cualquier sección nos referimos a los temas por medio de un número. Al referirnos al mismo tema en otra sección, citamos 
cuatro números: número de la unidad, número del capítulo o parte, número de sección y número del tema. Por ejemplo, “Ecuación 3” 
si aparece en la misma sección, y “Ecuación 4.3.1.3”, si está en otra sección. Los tres primeros números se citan en la parte superior 
de la página. Creemos que este es un sistema claro y fácil para consultar cualquier tema. Los indicadores serán una guía valiosa 
y le ahorrarán tiempo. 

En la margen derecha del texto hemos tratado de indicar la clase del material. Los términos usados son tan claros que no necesitan 
explicación. Por ejemplo, “Discusión” significa que una vez que haya entendido el tema, no necesita de nuevo recurrir a esta parte 
del texto. Estos INDICADORES tienen por objeto ayudarlo a localizar rápidamente las partes importantes cuando esté repasando 
o tenga que omitir parte del texto por falta de tiempo. 

También hemos tratado de indicar la importancia de los temas por el sistema de estrellas. 

*** Esto es muy importante. Usted debe entenderlo claramente. En algunos casos lo usará con tanta frecuencia, que se familiari- 
zará con ello sin necesidad de memorizarlo deliberadamente. El punto principal de este curso es no aprenderse los datos de me- 
moria; nosotros creemos que es más importante que entienda el trabajo y sea capaz de aplicarlo en la solución de los problemas. 
** Este es de menor importancia, pero tiene material al cual nos referiremos más adelante. 

* Esto generalmente no lo guía a trabajos adicionales en este curso y lo puede omitir si tiene poco tiempo. 

M Indica el final de cada ejemplo, ejercicio o solución. 
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Objetivos 


El objetivo principal de esta unidad:es introducir las ecuaciones di- 
ferenciales de primer orden y discutir sus soluciones y algunos méto- 
dos para resolverlas. 


Después de haber trabajado esta unidad debe estar en capacidad de: 


(i) explicar el significado de los siguientes términos: 
conjunto solución de una ecuación diferencial, 
familia de curvas solución, 
método de separación de variables, 
método del factor de integración, 
método de Euler para resolver ecuaciones diferenciales de pri- 
mer orden; 

(ii) hacer un diagrama de una familia de curvas solución de una 
ecuación diferencial dada de primer orden; 

(iii) resolver una ecuación diferencial dada de primer orden usando 
el método de separación de variables o el método del factor de 
integración (en los casos apropiados); 

(iv) una vez completada la parte (iii), encontrar una solución particu- 
lar con una condición inicial dada; 

(v) aplicar el método de Euler, dadas una ecuación diferacñal de 
primer orden, una condición inicial y una longitud de paso es- 
pecificada. 
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Glosario 


CONDICION INICIAL 
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PRIMER ORDEN 


CONJUNTO SOLUCION 
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DIFERENCIAL 


CURVA SOLUCION 
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LONGITUD DE PASO 
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Una condición inicial de una ecuación dife- 
rencial de primer orden es una condición que 
selecciona un elemento particular del con- 
junto solución. Geométricamente representa 
un punto dado, a través del cual pasa la grá- 
fica de ese elemento particular. 


El conjunto solución de una ecuación dife- 
rencial es el conjunto de todas las funciones 
reales que verifican la ecuación diferencial. 


Una curva solución de una ecuación diferen- 
cial es la gráfica de un elemento del conjunto 
solución. 


Una ecuación diferencial es una ecuación 
que relaciona una función real desconocida 
(que queremos encontrar) con sus funciones 
derivadas y, posiblemente, con otras funcio- 
nes conocidas. Tal ecuación se puede expre- 
sar en términos de las imágenes (variables). 
(Véase en el texto la discusión detallada.) 


La familia de las curvas solución de una 
ecuación diferencial es una representación 
gráfica del conjunto solución. 


Una función de dos variables reales es una 
función cuyo dominio es un subconjunto de 
R X R, y su dominio es R. 


Una función de una variable real es una fun- 
ción cuyo dominio es un subconjunto de R, 
y cuyo codominio es R. 


La longitud de paso es la longitud del inter- 
valo usado para determinar el punto siguien- 
te de una curva solución estimada de una 
ecuación diferencial de primer orden. 


El método de Euler es un método numérico de 
solución de una ecuación diferencial. (Véa- 
se en el texto la discusión detallada.) 


Un método de fórmulas es un método de so- 
lución de una ecuación diferencial que pro- 
duce una fórmula para la solución. 


El método de separación de variables es un 
método particular de solución de una ecua- 
ción diferencial de primer orden. (Véase en 
el texto la discusión detallada.) 


El método del factor de integración es un mé- 
todo especial de solución de una ecuación di- 
ferencial de primer orden. (Véase en el texto 
la discusión detallada.) 
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ORDEN DE UNA 
ECUACION DIFERENCIAL 


SOLUCION DE UNA 
ECUACION DIFERENCIAL 


VARIABLE 
INDEPENDIENTE 


El orden de una ecuación diferencial es el 
orden de la mayor función derivada que apa- 
rece en la ecuación. 


Una solución de una ecuación diferencial 
es un elemento del conjunto solución. 


Si Q:t-—>q, entonces t, la variable del do- 
minio de Q, se llama variable independiente, 
y q se llama variable dependiente. 
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Bibliografía 


La mayoría de los libros de textos son más profundos que este texto. 
Por esta razón hemos escogido solo 2 libros para darle alguna idea del 
número de métodos clásicos útiles, y para dar un punto de vista dife- 
rente a algunas partes del material presentado en este texto. Ellos son: 


H. Betz, P. B. Burcham y G. M. Ewing, Differential Equations with 
Applications (Harper and Row, 1964). 
Son importantes las partes de los capítulos 1-4. 


W. Kaplan, Elements of Ordinary Differential Equations (Addison- 
Wesley, 1964). 
Son importantes las partes de los capítulos 1-3. 


Para el enfoque numérico sugerimos: 


N. Macon, Numerical Analysis (John Wiley, 1963). 

El capítulo 10 da una indicación del enfoque más formal de los méto- 
dos numéricos y su exactitud. Usted encontrará probablemente esta 
última referencia más bien difícil en esta etapa. 


Nota general 


En esta unidad usaremos cualesquiera de las notaciones del cálculo 
que hemos mencionado previamente en este curso, según nos parez- 
can, en cada caso, apropiadas o convenientes: Por ejemplo, si Q es 
una función real tal que 


O:tH—=>0(1) (te R) 


y escribimos q = Q(t), entonces denotaremos la función derivada por 
Q', cuando la “notación de función” nos convenga, o expresaremos la 
derivada por 
dq 

AC 
usando la “notación de imagen de la función” o la notación de Leib- 
niz. Ya debe estar familiarizado con la notación de función puesto que 
la hemos usado a través de todo el curso y, también, con la diferencia 
entre la función Q y las variables t y q, donde t es un elemento gene- 
ral del dominio de Q, y q es un elemento general del conjunto imagen. 
Nuestra intención, al usar estas distintas notaciones, es capacitarlo 
no solamente para que lea este curso, sino también para que pueda 
leer cualquier otro texto que encuentre ahora o en sus estudios futuros. 
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24.1 INTRODUCCION 


24.1.0 ¿Qué es una ecuación diferencial ? 


En la Unidad 12, Diferenciación I, vimos que la derivada de una fun- 
ción real Q representa la razón de cambio del valor de la imagen en 
Q, o la pendiente de la (recta tangente a la) gráfica de Q; esto es 
+ Q'(t), (+ E R), es la razón de cambio de Q(t) en t. Dada la función Q, 

podemos encontrar Q' por medio de las reglas corrientes. Por ejemplo, 
si 

Oui —>e! = 21? (te R), 
entonces 

O':t— e! — 61? (tE R). 
Este es un método directo. 


El problema es mucho más difícil cuando nos dan la función derivada 
Q' y nos piden que encontremos Q. Este problema es la integración, 
que ya hemos discutido, y para la cual hemos desarrollado algunas 
herramientas (Unidad 9 y Unidad 13, Integración 1 y 11). Esencial- 
mente, hay dos puntos principales. En primer lugar, el operador deriva- 
ción D es una función, pero la operación inversa de D no determina 
una función. Para vencer esta dificultad introdujimos una función 
constante y obtuvimos las imágenes de Q' producidas por la aplicación 
inversa: 


cualquier primitiva + una función constante 
Por ejemplo, si nos dan 

O >+21 (tE R) 
entonces 

0i— + (1ER). 
donde c es un número real, y para cada escogencia de c obtenemos una 
imagen de (? en la aplicación integración. El segundo punto importan- 
te es que, con mucha frecuencia, no podemos encontrar una expresión 


sencilla para la integral (o, aun cuando se pudiera encontrar, el traba- 
jo requerido sería excesivo). 


Por ejemplo, si 


== (1ER*) 


entonces 


O: DE>2 


En esta unidad haremos más difícil aún el proceso de la integración; 
no porque ello ““sea formativo” sino por el amplio campo de aplicación. 
Estamos en los umbrales de la materia conocida por ecuaciones dife- 
renciales, sobre la cual existe abundante literatura y que, aún hoy, 
es objeto de profundas investigaciones. 


Una ecuación diferencial del tipo que consideramos en esta unidad 
es una función que comprende una función real desconocida y su fun- 
ción derivada, así como, posiblemente, otras funciones conocidas. Así, 
si f es una función real, 


PH =:0 
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Introducción 


k * * 


(donde O es la función x——0(x E R)) es una ecuación diferencial, 
como lo es 


f' + (xx?) x f = x—> sen x. 


La forma más sencilla de una ecuación diferencial es una ecuación de 
la forma 


f=g 


donde g es una función real conocida; su solución nos enfrenta preci- 
samente al problema de la integración que ya hemos discutido ante- 
riormente. En la Sección 24.1.2 discutiremos más ampliamente lo que 
entendemos por la solución de una ecuación diferencial pero, por aho- 
ra, vale anotar que la solución de una ecuación diferencial es una fun- 
ción (o un conjunto de funciones). Este hecho se diferencia de lo que 
sucedía en las ecuaciones que habíamos considerado hasta ahora, 
donde, generalmente, la solución era un número (o un conjunto de nú- 
meros). En general, una ecuación, en matemática, es una relación que 
está definida entre los elementos de un conjunto y que se puede con- 
siderar que define un subconjunto cuyos elementos son los mismos 
que pertenecen al conjunto solución de la relación. (Véase en la Uni- 
dad 19 la definición de relación.) En el lenguaje de la Sección 17.2.1, 
Unidad 17, Lógica II, una ecuación es una proposición abierta que se 
convierte en una proposición verdadera cuando la variable se sutitu- 
ye por un elemento del conjunto solución. 


Así, una ecuación diferencial es una relación definida en un conjunto 
de funciones reales y, por tanto, su conjunto solución está formado por 
funciones. Por ejemplo, un elemento del conjunto solución de la ecua- 
ción diferencial 


f+2f$=0 
es la función 
f:x=— -2exp(—2x) (x € R). 


Observe que una relación entre funciones se puede expresar en térmi- 
nos de las imágenes de las funciones. Por ejemplo, 


f+a=0 
es equivalente a 
FP) + 2(x) = 0, 


para todo x del dominio de f'. Nos referiremos también a la “forma ima- 
gen” de una ecuación diferencial como a una ecuación diferencial; en 
este caso, con frecuencia es más conveniente usar la notación de Leib- 
niz para las derivadas. 


Las ecuaciones diferenciales aparecen directamente de muchas leyes 
físicas básicas y, por consiguiente, son fundamentales para el estu- 
dio de gran parte de las ciencias y las ingenierías. Por ejemplo, una de 
las leyes de Newton establece que 


fuerza aplicada = masa X aceleración 


La aceleración es la razón de cambio de la velocidad, y la velocidad 
depende del tiempo. Si, por ejemplo, la fuerza aplicada depende de 
la velocidad, como sucede en el motor de un automóvil, entonces pode- 
mos escribir la ecuación anterior en términos de una función “posi- 
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ción” y de su primera y segunda derivadas, con lo cual obtendremos 
una ecuación diferencial. Las ecuaciones diferenciales no aparecen 
únicamente en esta clase de problemas: tienen otras muchas y diver- 
sas fuentes como la arquitectura, la biología y la economía. En esta 
unidad resolveremos solamente una forma sencilla de ecuación dife- 
rencial que contiene una función real Q, su función derivada Q”' y 
ciertas funciones conocidas. 


Nuestro plan consiste en.concentrarnos principalmente en unas cuan- 
tas ecuaciones diferenciales que estén bien relacionadas. En la Sec- 
ción 24.1.1 estableceremos el modelo; en otras palabras, desarrollare- 
mos algunas ecuaciones diferenciales tomadas de situaciones físicas 
aceptables. En la Sección 24.1.2 discutiremos algunas ideas básicas 
que se refieren a la solución de las ecuaciones diferenciales. Después 
investigaremos las soluciones gráficas de las ecuaciones diferenciales 
porque un esquema gráfico nos dará una buena idea del comporta- 
miento de las soluciones. Las soluciones gráficas nos permiten tam- 
bién apreciar el enfoque numérico que introduciremos en la Sección 
24.2.4. En las Secciones 24.2.2 y 3 introducimos lo que podríamos lla- 
mar métodos de fórmulas. El objetivo de un método de fórmula consis- 
te en reordenar la ecuación diferencial de modo que podamos usar 
nuestra poderosa herramienta: la integración. 


Este es el modo de enfocar la solución de un problema, según aconse- 
ja Polya. Nos preguntamos: “¿Hay algún problema análogo a este que 
ya hayamos resuelto?” La respuesta será un “sí” si comenzamos por 
hacer una pequeña reordenación. 


24.1.1 El problema del crecimiento 
de población 


En esta sección consideraremos el problema del crecimiento de pobla- 
ción que ya habíamos discutido en la Sección 7.4 de la Unidad 7, Su- 
cesiones y límites I, y examinaremos modelos matemáticos mucho 
más elaborados. 


Supongamos que q = Q(t) representa (en un tiempo cualquiera t E R) 
la población (en millares) de una especie particular, que considerare- 
mos que son seres humanos, pero que igualmente podrían ser virus, 
langostas, peces o pájaros. 


MB 24.1.0/24.1.1 


24.1.1 


Discusión 
h £*x 


Puesto que q se mide en millares, solamente podrá tomar ciertos va- 
lores racionales correspondientes a números enteros de seres huma- 
nos, digamos 285.632. Esto nos brinda la primera modificación en el 
proceso de formular nuestro modelo matemático. Para hablar de creci- 
miento o de razón de cambio, quisiéramos podernos referir a la fun- 
ción derivada, pero la función derivada no se ha definido para el ca- 
so de funciones cuyo codominio es este subconjunto de racionales. La 
razón es que el concepto de límite, que usamos en la definición de fun- 
ción derivada, carece de sentido en ese caso. Esto parece ser una ba- 
rrera infranqueable pero, para saltarla, tomamos como codominio de 
Q el conjunto R*, de modo que podamos construir nuestro modelo ma- 
temático. Esto es, suponemos que la variable q puede tomar cualquier 
valor real positivo, de manera que podamos utilizar las poderosas he- 
rramientas del cálculo. 


Para determinar cómo cambia la población con el tiempo necesitamos 
introducir funciones para representar el número (en millares) de naci- 
mientos y de muertes anuales. Las denotaremos por B y M, respecti- 
vamente. El número de nacimientos (y de muertes) depende general- 
mente del tamaño de la población, Q(t), que, a su vez, depende del 
momento en que se considere. Así, la función compuesta B + Q nos 
dirá de qué manera el número de nacimientos anuales depende del 
tiempo. 

El dominio de cada una de las dos funciones B y M es el subconjun- 
to de números reales que representa la población en millares, y el co- 
dominio es un subconjunto de los reales que representa el número de 
nacimientos (o muertes) anuales. Mediremos t en años. La razón de 
crecimiento viene dada, entonces, por 


DQ =B-Q-M-Q 


Q'(1) = B(Q(t) — M(Q(1)) = Blg) — M(q) 


Lo más sencillo que podemos suponer es que el número de nacimien- 
tos y de muertes que ocurren anualmente son constantes; es decir, 
B:ig=—>b9  (qeR”), 
M:q3-—mMo (qeR?), 
donde bo, y my son números conocidos. 
Así, obtenemos 
01) = by — My = ky. por ejemplo 
No es nuestro propósito resolver ecuaciones diferenciales en esta sec- 
ción, aunque ésta tiene una solución obvia. Lo que queremos ilustrar 
aquí es cómo aparecen las ecuaciones diferenciales. Ahora bien, es 
irreal suponer que el número de nacimientos anuales permanecerá 
siempre constante; por ejemplo, si la población se duplica, es válido 
suponer que el número de nacimientos también se duplicará. Para 


trasladar esta idea a términos matemáticos necesitamos que la fun- 
ción B sea de la forma 


Bs >b:g (qeR?), 


donde b; es un número positivo. Esto es, el número de nacimientos 
es proporcional al número de sujetos que forman la población. 
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Del mismo modo, es razonable suponer que el número de muertes es 
proporcional a la población; es decir, 


M:q——m:q (qeR?), 


donde m, es un número positivo. La ecuación diferencial es, ahora, 


Q'(1) = b,0(1) — m,Q(t) = k,0(t) 


—= kg, 


donde k; = bi; — mi. 


Ahora refinaremos nuestro modelo. Como la población crece, la canti- 
dad de alimentos disponibles comenzará a disminuir. Con una mayor 
posibilidad de enfermedades y accidentes que provocan un mayor nú- 
mero de muertes, debemos esperar que la razón de muertes ha de cre- 
cer con el crecimiento de población. 


Esto lo podríamos, quizá, representar matemáticamente por 
M:q——>m q + mq? (qeR*) 


donde m2 es un número positivo. Por ejemplo, si tomamos m, = 0,02 
y ma = 0,004, tendremos 


M:q-—>0,02q + 0,004q? (q E RF) 


la gráfica de M se muestra en la figura: 


4 M(a) 


0,6 7 


0,5 y 


qr»0,02q + 0,004 q2 


Esto implicaría que, con una población de 1000 (esto es, con q = 1), 
la razón de muertes debería ser de 24 por mil por año, mientras que, si 
la misma población aumentara a 10.000 (esto es, si q = 10), la ra- 


0,2 + 0,4 


zón de muertes aumentaría a = 0,06, es decir, 60 por mil por 


año. La ecuación diferencial sería, ahora, 


0'(1) = k,Q(1) — ma(Qlo))? 


donde k; = b; — m,. 


De ahora en adelante escribiremos la “forma imagen” de una ecuación 
diferencial en términos de Q(t) o q, en vez de continuar usando ambas 
versiones. La primera forma, explícitamente en términos de las imá- 
genes, es útil porque nos dice que como solución tenemos una función 
en vez de un número. La segunda forma, con la forma de Leibniz, es 
conveniente porque es más concisa. Usaremos una u otra forma de 
acuerdo con lo que nos parezca que es más apropiado al contexto. 


Es interesante observar que una ecuación diferencial nos brinda al- 
guna información aun cuando no la hayamos resuelto. Podemos es- 
cribir nuestra ecuación diferencial en la forma 


dq dl k; 

> 24 m, q 

dr , 2 
po sería negativa y la población 
decrecería hasta que el modelo resultara inaplicable o la población lle- 
gue a cero) vemos que, si 


Suponiendo que k, > 0 (si no lo fuera, 


As positiva, lo cual implica un crecimiento de población. (Recuerde 


dt 
6 
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: . 4 y k 
que q = Q(t), es decir, q es un número.) Si q a entonces la po- 
E a e k z Z e 
blación es estática, y si q >—, la derivada es negativa y la población 
m, 


k 
decrece. Esto nos dice que — es la única población estable. 


2 
Esto nos lleva a la última refinación que pretendemos hacer aquí. Es- 
ta población estable que se puede mantener en esa forma, está llama- 
da a crecer con el tiempo, debido a las mejorías que se logran en la 
medicina, en la tecnología agrícola, etc. Supongamos que esta pobla- 
E , : E UÍn k 
ción estable crece linealmente con el tiempo y que sustituimos — por 


mM) 
(Rat + k3), donde k2 y kz son números positivos. La ecuación di- 


ferencial se convierte en 


dq 
5 de m2Q (Rat + kz — q) 

Observe que ahora el miembro de la derecha comprende las dos varia- 
bles q y t. A estas variables se les da generalmente nombres especia- 
les: t, la variable del dominio de la función Q que queremos determi.- 
nar se llama Variable independiente, q es la Variable dependiente 
porque es la variable del codominio de la función Q que queremos de- 
terminar. 


Observe también que hemos supuesto algunas cosas sin que tuviéra- 
mos ninguna justificación verdadera. Lo que tenemos es un modelo 
tentativo que, sometido al análisis, puede ser que se adapte a los he- 
chos o no. En una situación real necesitaríamos comprobar sus conse- 
cuencias valiéndonos de datos conocidos, antes de lanzarnos a sacar 
predicciones de él. 

Todas las ecuaciones diferenciales que tratamos en esta introduc- 
ción son ecuaciones de primer grado. El 0rden lo determina la mayor 
derivada presente. En otras palabras, si una ecuación diferencial con- 


In 


E d , ] ¿ 
tiene = O Q”(t) y no contiene derivadas mayores, decimos que es 


de orden n. Por ejemplo, la ecuación 


dig va 
(5 E 


es de orden 3, o de tercer orden. 


24.1.2 Ideas básicas acerca de las soluciones de 
las ecuaciones diferenciales 


Consideremos el problema familiar de resolver la ecuación cuadrática 
2-St+44=0  (teR). 

El conjunto solución de esta ecuación, a saber 
(t:1? — St + 4=0,teR), 

contiene dos elementos y es 


(4,1). 
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Definición 1 
h $ $ 


Definición 2 
h $x e 


24.1.2 


Discusión 
kk 


El conjunto solución, en R, de una ecuación cuadrática puede conte- 
ner dos elementos, un elemento (cuando las dos raíces son iguales) o 
ningún elemento. 


Ejercicio 1 


Escriba ejemplos de tres ecuaciones cuadráticas tales que los conjun- 
tos soluciones contengan, respectivamente, dos, uno y ningún ele- 
mentos reales. un 


Ejercicio 2 
Describa en forma explícita los elementos del conjunto solución de la 
ecuación 

sent=0 (te R) 


¿Puede contar cuántos elementos tiene ese conjunto? E 


El conjunto solución de una ecuación diferencial se puede representar 
de una manera semejante al de una ecuación ordinaria. Considere la 
ecuación diferencial 


Q()=k,  (teR) 


010) — ko =0 


El conjunto de todas las funciones que verifican esta ecuación dife- 
rencial se puede denotar por 


(0:00) -k,=0  (teR)). 


En este caso cada solución de la ecuación diferencial es sencillamen- 
te una función primitiva (integral indefinida) de tt——ko, esto es 


tt kot+<c (te R); 


c es una constante de integración (véase la Unidad 13, Integración 
11). Así, otra forma explícita del conjunto solución es 


[0:0 = tt kot+< (teR),ceR), 


que con frecuencia se expresa más concisamente en términos de imá- 
genes como 


(Q:0(t) = kot + c (teR),ceR). 


Este conjunto también se puede escribir como una relación entre va- 
riables: 


(0:q =kot+c<c (teR),ceR). 


De nuevo podemos observar aquí la importante diferencia que existe 
entre las soluciones de las ecuaciones ordinarias definidas en conjun- 
tos de números y las soluciones de las ecuaciones diferenciales defi- 
nidas en conjuntos de funciones. En el primer caso el conjunto solu- 
ción tiene, como elementos, números reales mientras que en el se- 
gundo caso tiene funciones reales. Visualmente podemos observarla 
si consideramos la representación gráfica del conjunto solución de la 
ecuación 


sent=0 (te R), 
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Ejercicio 1 
(2 minutos) 


Ejercicio 2 
(2 minutos) 


Discusión 
* * 
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y la comparamos con la representación gráfica del conjunto solución 
de la ecuación diferencial 


O()-1=0  (teR) 

Por el ejercicio 2 sabemos que el conjunto solución de 
sent=0 (tE R) 

es el conjunto de números 
Ínr:ne Z). 


Esto se puede representar gráficamente como la intersección de la rec- 
ta que es representación de 


t—>0 
y la curva que es representación de 


t>Sent. 


Las coordenadas de los puntos de intersección son de la forma (nr, 0), 
DAA 


sen t 


trHesen t 


El conjunto solución de la ecuación diferencial 
Q()-1=0  . (teR) 


se puede representar gráficamente por medio de la familia de curvas, 
o conjunto de curvas solución, que viene dada por la representación 
gráfica de las funciones 


id a (teR),ceR. 
La causa de que exista una familia proviene de todos los diferentes 


valores que puede tomar la constante de integración c. Ilustramos los 
casos donde c = 2, 1,0, —0,6, —1, —2. 


Todo el plano queda cubierto por curvas que, en este caso, son líneas 
rectas. Por un punto dado cualquiera pasa una curva de la familia, y 
solamente una. Seleccionar un punto del plano equivale a seleccionar 
una curva particular, lo cual es equivalente a seleccionar una constan- 
te particular de integración. Por ejemplo, consideremos una población (continúa en la página 10) 


9 


Solución 1 


Por ejemplo, 


(i) tz — 3t + 2 = 0 cuyo conjunto solución es (2, 14. 
(ii) tz — 2t + 1 = 0 cuyo conjunto solución es fi 
(li) ta — 2 + 2 = 0 cuyo conjunto solución no contiene números 
reales. ul 
Solución 2 


El conjunto solución se puede escribir en la forma 
(nm:ne Z), 


porque el seno de cualquier múltiplo entero de r es cero. Hay tantas 
soluciones como números enteros positivos. Las soluciones se pueden 
poner en correspondencia biunívoca (uno a uno) con el conjunto de 
todos los enteros positivos como sigue: 


y, generalmente, 
kn ——> 2k para k > 0, 
kx ——> —2k + 1 para k <0. | 


(viene de la página 9) 

de peces que comienza, cuanto t = 0, con 2000 miembros y que crece 
con una razón constante de 1000 por año (este número representa la 
diferencia entre el número de nacimientos y el número de muertes). 


Q(t) 


Entonces, usando unidades de 1000, tenemos un punto inicial, (0, 2), 
(señalado con una X en el diagrama), y la curva solución apropiada 
(señalada con un trazo de color en el diagrama) es la gráfica de 


0 (te R) 
puesto que nuestra solución general (i.e. Q:t——>t + c) implica que 
Q(0)=0+c, 
1.6, 
2=0+c. 
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Solución 1 


Solución 2 
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En las ecuaciones diferenciales que se obtienen como modelos mate- 
máticos de situaciones reales, generalmente conocemos un valor ini- 
cial, tal como vimos en el ejemplo de los peces. Aquí, si el modelo es 
satisfactorio, la curva solución particular que pasa por el punto 


(tiempo inicial, población inicial) 


se puede utilizar para determinar la población futura. Este proceso de 
escoger una curva particular por medio de la condición inicial conoci- 
da es una característica importante del modelo matemático que se 
ha construido para esta situación. 


Ejercicio 3 Ejercicio 3 
” ' e E (3 minutos) 
Represente gráficamente un subconjunto típico de la familia de curvas 


que ilustran el conjunto solución de la ecuación diferencial 
Qt) - 2t =0 (te R). 


Indique la curva solución particular que verifica la condición inicial 


Q(1) = 2. ¿Cuál es, en ese caso, el valor numérico de la constante 
de integración? E 
Ejercicio 4 Ejercicio 4 
; (2 minutos) 
Escriba 
Q'(t) = 2 
en una forma que comprenda 
(i) funciones solamente. 
(ii) variables solamente. n 
Las ecuaciones diferenciales, como las ecuaciones ordinarias, pueden Discusión 
ser tales, en algunos casos, que los correspondientes conjuntos solu- sd 
ciones sean vacíos. Así, 
LAUIOY +00 +1=0 (te R)) 
es vacio, del mismo modo que 
(tu? -21+2=0 (te R)) 
es vacio. Ya hemos visto que el número de soluciones de Q'(t) — 2t =0 
no es finito. También una ecuación diferencial puede tener un núme- 
ro finito de funciones en su conjunto solución. Así, 
UFO? + (0) =0 (te R) 
tiene un elemento, 
Jet 0 (tE R), 
del mismo modo que 
(t:1?-21+1=0 (te R)) 
tiene solamente un elemento. 
Resumen 


Resumen 


* 
Las soluciones de las ecuaciones diferenciales son funciones. Las 
gráficas de estas soluciones forman una familia de curvas solución. 
En nuestro ejemplo usamos una condición inicial (un punto del plano) 
para escoger una función particular (una curva solución particular). 
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Solución 3 


tht2+2 


lA 
trot2+1 


El conjunto solución es 
10:01) =1*+c (teR),ceR). 
Utilizando la condición inicial 
Q(1l)=2=1+c 
obtenemos 
c=1 
de donde la curva solución particular corresponde a 
Q(t) = 1? +1. 


La familia de curvas es, en este caso, una familia de parábolas. Todas 
las curvas parecen semejantes en el mismo sentido que una familia 
biológica puede tener rasgos semejantes. Esta es una razón por la cual 
es conveniente usar la palabra familia para describir esta colección 
particular de curvas. [E] 


Solución 4 
(1) 0” =1f=21 (te R) 
o  DQ=tr—>2t (te R) 
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Solución 3 


Solución 4 


24.2 ALGUNOS METODOS PARA RESOLVER 
ECUACIONES DIFERENCIALES 


24.2.0 Introducción 


En esta sección analizamos tres modos de encontrar las soluciones de 
las ecuaciones diferenciales de primer orden. Cada uno tiene sus ven- 
tajas y sus desventajas, que iremos señalando en la discusión. 


En la Sección 24.2.1 discutimos un método gráfico de solución y da- 
mos un ejemplo. 


En las Secciones 24.2.2 y 3 seleccionamos, del extenso repertorio de 
métodos de fórmulas, dos de ellos que dan soluciones exactas de ecua- 
ciones diferenciales de primer orden de ciertos tipos. (Cuando deci- 
mos que las soluciones son exactas, queremos decir que podemos es- 
pecificar con precisión el conjunto de funciones que forma el conjunto 
solución, en vez de dar una tabla de soluciones como se suele hacer 
en un método numérico.) La utilidad de los métodos de fórmulas con- 
siste en que producen una solución más “general” que el método grá- 
fico, del mismo modo que las soluciones de la ecuación cuadrática 


ax?+bx+c=0 
se expresan de una manera más general por medio de la fórmula 


=b+./b? — 4ac 


2a ] 
que por medio de algunas condiciones, tales como las de hacer a = 1, 
b= —2yc = 1, para concluir que, entonces, x = 1. La solución ge- 


neral se puede evaluar para muchos valores numéricos de a, b y c. 


¿Por qué discutimos solamente dos métodos de fórmulas? La razón es 
que las dos que hemos escogido son fáciles de aprender, y son repre- 
sentativos de los métodos exactos. No es nuestra intención hacer una 
lista de todas las posibles fórmulas. En vez de eso, esperamos que ob- 
tenga una idea general de lo que significa buscar soluciones exactas 
de ecuaciones diferenciales, y que se capacite para poder aprender 
otros métodos cuando sean necesarios. 


Todos los métodos exactos descansan necesariamente en los métodos 
de integración y derivación que discutimos en las Unidades 12 y 13. 
Los mencionaremos explícitamente cuando se necesiten. 


Por último, en la Sección 24.2.4 damos un ejemplo de un método numé- 
rico de solución. 


24.2.1 Métodos gráficos 


Pudimos escribir directamente la solución de la ecuación 

sent=0 
porque sabíamos que 

sen (nx) = 0 (ne Z), 
y que, además, ningún otro valor de t del dominio R de la función se- 
no está aplicado a cero. Ahora bien, supongamos que tenemos que re- 
solver una ecuación de la forma f(t) = O que resulte: menos fácil; por 
ejemplo, la ecuación que teníamos en “el problema de la tortilla”, pá- 
gina 34 de la Unidad 2, Errores y exactitud: 


2x1 
t— t=-—=0, 
sen 3 0 
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24.2 


24.2.0 


Introducción 
* $ 


24.2.1 


Discusión 
* * 


Allí utilizamos métodos gráficos. Calculamos las imágenes de la fun- 
ción 


NH =sent E (te R) 


en algunos puntos especiales del dominio, los unimos por medio de un 
trazo suave (suponiendo que la función es continua); el valor de t don- 
de la curva trazada cortó el eje x era la solución aproximada que bus- 
cábamos. 


pp 


(3,14,1,05) 


(2,09,-0,87) 


(1,05,-1,91) 


0,-2,09) 


Gráfica aproximada de tre t-sent- 271 


Para mejorar la exactitud, ampliamos la región de la curva en los alre- 
dedores de la solución aproximada calculando nuevas imágenes de la 
función y trazando esa parte de la gráfica más precisamente. 


Podemos adoptar un procedimiento semejante para buscar la solución 
de las ecuaciones diferenciales de primer orden; este procedimiento 
será muy útil en aquellos casos que no podamos resolver por medio de 
los métodos comunes que discutimos en el texto. La última ecuación 
del ejemplo del crecimiento de población de la Sección 24.1.1: 


dq 
FA maq(x2t + k3— 9) 
t 
es uno de esos casos. Usaremos, sin embargo, la ecuación previa 
dq 
— ="k,q — m>q? 
dt 19 24 


para ilustrar el método gráfico. Es obvio que para lograr una solución 
gráfica, debemos usar números en vez de k, y mo». Esto ilustra la li- 
mitación de los procedimientos numéricos o gráficos; no solamente te- 
nemos que especializar el problema aplicándolo a casos particulares 
sino que de ninguna manera obtendremos una expresión general de 
la solución. Supongamos que escogemos k, = 2 y m2 = 1. La ecua- 
ción diferencial es, entonces, 

dq 3 

q 974 
y queremos encontrar la función, o funciones, Q, tales que q = QÍt), 
para verificar esa ecuación. Nos contentaremos con tener tal función, 
o funciones, especificadas gráficamente. A cada una de las correspon- 
dientes curvas la llamaremos curva solución. 
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E a AA dl 
Por la Unidad 12, Diferenciación I, sabemos que = = (Q' (t) repre- 
senta la pendiente de la gráfica de la función Q; esto es, representa la 
dirección hacia donde apunta la curva solución en un punto cualquie- 
ra (t, q). Para nuestra ecuación diferencial podemos calcular la pen- 
diente de una curva solución en cualquier punto; siempre viene dada 
por 


2q — q?. 
Así, en el punto (0, 1), que representa una población de 1000 en el 
tiempo t = O, la pendiente de la curva solución particular es 


2x1-1*=1 


Puesto que 2q — q? no depende explícitamente de +, la pendiente 
en cualquier punto (t, 1) es siempre 1. En términos del crecimiento de 
la población que representa la ecuación diferencial, esto significa, 
como esperaríamos, que la razón de crecimiento depende explícita- 
mente del tamaño de la población y no del tiempo en que se esté mi- 
diendo. En términos de las curvas soluciones, esto significa que, para 
q > 0, tenemos una curva solución que pasa por todos los puntos (t, q), 
y que, para cada q, las rectas tangentes a las curvas soluciones en 
(t, q) son paralelas entre sí, cualquiera que sea el valor de t. 


Señalemos las pendientes en un conjunto de puntos, mediante una 
gráfica apropiada. Recordemos que la pendiente es la razón 


distancia vertical 
distancia horizontal 


Por consiguiente, fijemos una distancia horizontal de, digamos, cua- 
tro cuadrados del papel cuadriculado. Después, tracemos la flecha que, 
en un punto particular, nos indica la dirección de la curva solución. 
Por ejemplo, para el punto (1, 0,5) necesitamos trazar una flecha que 
represente una pendiente de 0,75. El diagrama indica cómo se hace 
esto. 


En la tabla que ofrecemos a continuación damos los valores numéricos 
verdaderos de las pendientes en algunos otros puntos seleccionados. 
Como hemos observado, el valor de la pendiente no depende del valor 
de t; en consecuencia, todas las columnas son iguales. No obstante, 
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hemos construido una tabla por dos motivos: primero, para que tenga 
una idea de lo que ocurre en un caso más general y, segundo, para ob- 
tener un registro de datos que es semejante a la gráfica. 


EEES 


Si trazamos las flechas correspondientes, obtenemos el diagrama si- 
guiente: 


EEES 
es A 

mu o 
eee 
a 
E HE 2 2 ms a : 


mE 
la| 


+++++ 


103 
an 
| 
HU 


Hay varias observaciones acerca de la tabla y del diagrama: 


(i) Puesto que todos los números que aparecen en cada fila son 
iguales, tenemos que dondequiera que las curvas soluciones 
corten una recta horizontal dada deben estar apuntando en la 
misma dirección. 

(ii) Normalmente nos interesa únicamente la parte del plano donde 
q >0yt > 0 (es decir, el primer cuadrante) puesto que esta es 
la región del plano donde tiene sentido este problema. 
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(i1i) Podríamos ir calculando la pendiente de más y más puntos del 
primer cuadrante y hacer el conjunto de flechas tan numeroso 
y denso como quisiéramos (limitados solamente por el ancho y 
el largo de los trazos del lápiz). 

(iv) Del mismo modo que podemos trazar una curva valiéndonos de 
unos cuantos puntos, podemos, también, trazar una curva, O 
curvas, usando las flechas como guías. El siguiente diagrama 
muestra algunas curvas soluciones típicas. Todo lo que tene- 
mos que hacer para trazar una curva solución es comenzar en 
un punto dado y asegurar que la dirección de la curva se vaya 
manteniendo paralela, en cuanto sea posible, a las flechas que 
estén cerca de ella. 


(v) Si queremos mejorar la exactitud, debemos usar más flechas, 
del mismo modo que usaríamos más puntos para trazar una 
gráfica corriente. 

(vi) Lo que hemos obtenido es una aproximación de una selección 
de curvas de la familia de curvas que representa la solución 
de la ecuación diferencial. Hemos obtenido, aunque laboriosa- 
mente, un esquema de la “forma” de la familia de las curvas so- 
lución. Si se observa de lejos este conjunto de flechas, captamos 
la impresión del “flujo” de las curvas. 

(vii) Para nuestra ecuación diferencial particular, tenemos que qa=0 
es una solución trivial; si la población inicial es cero, el razona- 
miento matemático, como otro cualquiera, nos dice que la po- 
blación será siempre igual a cero. Cuando q > 0, todas las cur- 
vas soluciones parecen aproximarse al mismo valor numérico, 2, 
cuando t se hace grande. Así, 2000 es el número de la población 
estable. Si inicialmente hubiera más de 2000, se reduciría a es- 
te número; si menos, aumentaría hasta ese número. Como men- 
cionamos en la última sección, podríamos haber deducido esta 
información particular sin necesidad de resolver la ecuación di- 
ferencial. En efecto, la derivada, Q'(t), es cero cuando 


29 — q? =0, 


esto es, cuando q = O (y hemos eliminado este caso) y cuando 
q = 2. El valor q = 2 nunca se alcanza en nuestro modelo, ex- 
cepto en el caso particular en que la población comience con 
2000. 

(viii) Para seleccionar la curva solución particular que representa un 
caso particular de interés, necesitamos más información; por 
ejemplo, una condición inicial. Supongamos que comenzamos 
con 750 de población; esto es, q = 0,75, y comenzamos a medir 
el tiempo desde ese instante, de manera que q = 0,75 cuando 
t = 0. Comenzamos a trazar la curva desde ese punto y a pro- 
ducir la curva solución aproximada que, en el diagrama ante- 
rior, se indica con un trazo punteado en color. Después de 2 
años* vemos (según nuestra curva) que casi se ha alcanzado 
la población de 2000. Estimamos que, al cabo de un año, debe 
haber una población de 1700. 


El método gráfico sirve para todas las ecuaciones diferenciales de pri- 
mer orden que podamos poner en la forma 


d 
e J(.g), 


dt 
con cierta condición impuesta a f. Todo lo que hay que hacer es eva- 
luar las imágenes producidas por f en puntos apropiados del plano. 
Las condiciones que, impuestas a f, aseguran que existe una solución 
son muy complicadas para este Curso Básico y las trataremos en un 
curso posterior. 


El método gráfico puede ser muy útil en los casos en los cuales no se 
pueda alcanzar la solución por otros medios. Puede ser inexacto (como 
cualquier trazado de curva de todo tipo) y, verdaderamente, es muy 
trabajoso pero tiene la ventaja de que ofrece una idea cualitativa de 
las formas de todas las curvas soluciones, esto es, del conjunto solu- 
ción como un todo. Con esta visión “panorámica” podemos escoger 
puntos particulares de interés y encontrar numéricamente las curvas 
soluciones particulares que pasan por ellos. Otra ventaja del método 
gráfico es que brinda una introducción a los métodos numéricos de so- 
lución de ecuaciones diferenciales. Esto lo veremos en la Sección 
24.2.4. 


Resumen 


En el procedimiento gráfico usamos la información geométrica que nos 
da la ecuación diferencial 


dq 
o h » 
a FM) 
a saber, las pendientes de las curvas solución en todos los puntos 


(t, q) del dominio de f, para trazar la familia de las curvas soluciones. 


24.2.2 Método 1 de fórmulas: separación de 
variables 


Este método es útil solamente cuando se pueden separar las variables 
de la ecuación. 


*¡Más vale tarde que nunca! 
¿Fue usted más rápido que el autor? 
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Resumen 
* * 


24.2.2 


Tema principal 


k hx $ 


Iustraremos el método con un ejemplo. Consideremos la ecuación 


0) =1Q(1)  (teR). 
Tomamos t y Q(t) = q como variables, y podemos separarlas reorde- 
nando la ecuación en la forma 

2 _, 

Ql:) 


supuesto que excluimos del dominio de Q cualquier t tal que Q(t) = 0. 
Escribiendo esto en forma de función (en vez de imagen) obtenemos 


IO e pt (te R). 


Q 


1 
Puesto que 0 Xx DQ y t——t son funciones iguales, podemos igualar 


sus primitivas, supuesto que escogemos correctamente las constantes 
de integración. Por consiguiente, tratamos de integrar ambos miem- 
bros. 


Podemos integrar el miembro de la derecha para obtener el conjunto 
2 


» os t 4 
de funciones primitivas de la forma > + c, C ER, pero la inte- 


gral del lado izquierdo puede no ser tan obvia*. En la Unidad 13, 
Integración II (Sección 13.2.4) obtuvimos la siguiente fórmula para la 
primitiva de una función compuesta: 


[or x Dk = [fs)os 


Comparándola con nuestra ecuación, vemos que necesitamos escoger 


1 
k=Qygo>ek e das. 


Entonces, según nuestra tabla de integrales, tenemos que 


fe=o—=t00+<, (teR*) 


[fe ok=t=>InQ()+c,  (QUER”), 


donde c, es algún número real. Observe cómo ha surgido aquí otra 
condición, Q(t) > O (porque el dominio de la función logaritmo es R*). 
Podríamos trabajar con Q(t) < 0 si fuera de algún interés, pero tendría- 
mos que hacerlo separadamente. Cuando usemos resultados normales, 
como hacemos aquí, debemos tener siempre cuidado de usarlos en las 
circunstancias apropiadas. 


Reuniendo las dos partes tenemos que 


? 
tt ——> In Q(t) + c, = pS +c 


*De hecho, estamos en un terreno familiar; véase el ejercicio 2, Sección 12.3.2 de la 
Unidad 12. 
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Esto especifíca el conjunto de funciones Q que pertenecen al conjunto 
solución de la ecuación original pero en una forma inconveniente: la 
variable dependiente q = Q(t) no está explícitamente expresada en 
términos de la variable independiente t. 


Así, pues, tratamos de hacer una reorganización. Primero usamos la 
forma más sencilla en términos de imágenes: 


p? 
Ing+c=3+c 
2 
1.€. 
p? 
Ing => +(c-—c;) 
2 
Puesto que c y c, son números reales cualesquiera, c — c, es tam- 
bién un número real. Podemos, por consiguiente, usar d E R para de- 
notar c — c, o, sencillamente, eliminar c,. (En general, no necesi- 
tamos introducir la constante de integración en ambos miembros de 
la ecuación puesto que ya sabemos que, si dos funciones son iguales, 
entonces sus primitivas se diferencian en una constante.) Por tanto 
ahora tenemos que 
2 


t 
Ing=>5+d : 


Recordando que ln es una función uno a uno y que exp es su inversa, 
tenemos que 


p? 
q = exp [5 a ¿] 


Ahora tenemos la variable dependiente q = Q(t) expresada explícita- 
mente en términos de la variable independiente t, y podríamos decir 
que tenemos una solución satisfactoria de nuestro problema. En rea- 
lidad podríamos arreglar las cosas un poco más si recordamos que 


*x 


expx=e 


Entonces, 
tE 2 2 
exp [5 +d| = e 1V+4 = e? x el 


Si hacemos A = e”, la solución general se convierte en 


p? 
q=4e"?0 Aexp (5) 
donde A E R*, puesto que exp d no puede ser negativa ni cero. 


La familia de las curvas soluciones se ilustran en la gráfica de la pá- 
gina siguiente. 


Hemos resuelto esta ecuación en pequeños pasos y vale la pena notar 
que hemos esparcido condiciones a lo largo del camino. Originalmente 
requeríamos que Q(t) % 0, y, más tarde, exigimos que Q(t) > 0, lo 
cual hace inútil la condición anterior. Por último, tenemos que 


p? 
Q(t) = Aexp (5) 


y, puesto que A E R*, la condición Q(t) > 0 se verifica automática- 
mente. Nótese también que las constantes de integración c y c,, que 
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. 
Gráficas de Q(1)=Aexp(-) 
para A=2,1 y 0,3 


originalmente aparecían añadidas al final de nuestras expresiones y 
que no tenían restricciones, se combinaron como resultado de nues- 
tras manipulaciones, de manera que la constante final, A, está condi- 
cionada (A > 0). Con un poco de práctica, el manejo formal de estas 
expresiones se hace decepcionantemente fácil pero, con la rutina, se 
suelen olvidar las consecuencias de lo que se está haciendo. 


Y eso, a pesar de haber tenido tanto cuidado. ¿Qué pasaría si no fué- 
ramos tan cuidadosos? Fácilmente se puede comprobar que 


p? 
Q(t) = A exp (5) 


donde A E R”, también pertenece al conjunto solución de la ecuación 
diferencial original. (Las curvas soluciones son semejantes a las des- 
critas en la última figura pero todas ellas están volteadas, por deba- 
jo del eje t.) Así, pues, el descuido trae sus consecuencias. Bien, eso 
puede suceder, igual que puede suceder en cualquier otro campo cien- 
tífico, pero debe ser un descuido disciplinado. En este caso, la disci- 
plina consiste en comprobar que nuestras manipulaciones formales, 
sin condiciones, no nos lleven a incluir en nuestro conjunto solu- 
ción algunas funciones que, de hecho, no son soluciones. 


Para generalizar el procedimiento que utiliza la separación de varia- 
bles necesitamos refundir la ecuación diferencial original en la forma 


(g20)x DO =h, 
o, en términos de variables (imágenes), 
dq 


y, después, esperamos poder integrar. Aquí está la razón del nombre 
del método; “separamos” todos los términos que contienen q y los es- 
cribimos en uno de los miembros de la ecuación para que podamos 
integrar. Esto significa que la ecuación diferencial original 


da 
q 169) 
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Ecuación (1) 


MB 24.2.2 


debe ser expresable en la forma 


a 
de gl) 
Ejercicio 1 Ejercicio 1 


A a > p (2 minutos) 
Existen algunas condiciones más refinadas que deben imponerse a g 


y h para que este método dé resultados positivos; no discutiremos en 
este texto esas condiciones, pero, ¿cuál es la restricción obvia que re- 
quiere g? ] 


Volviendo a la forma de la ecuación (1) tenemos que Tenia puncibal 


(g20) x DO =h ds 
Usando el resultado, que recordamos anteriormente, para integrar una 
función compuesta, obtenemos 


fe-oxm0=[1+0—0 


esto es, 


[fe -0- [1+0—0 


Así, si podemos obtener fe y fr de unas tablas apropiadas (quizá 


después de algunas manipulaciones), podremos resolver la ecuación 
diferencial. 


Resumen has 


Si + * 
D=-  (e0)+0,1eR) 
go0 Ñ o j 


entonces las soluciones verifican la ecuación 


[fe -0= fr+u—0 ceR. 


La forma de Leibniz de la regla de la separación de variables se sinte- 
tiza más adelante. 


Si 
dq ht) 


a 0), 
4 ne (g(q) 4 0) 


entonces las soluciones verifican la ecuación 


f soda = [nod + < 


Esta es una forma que, probablemente, encontrará en muchos libros de 
texto. Es posible que esta forma sea más fácil de recordar. 


22 


Ejercicio 2 


Encontrar la función Q que es solución de 


dq 
o O) R 
a t (teR) 


donde q = Q(t), y que satisface la condición inicial 
Q(0) = 3 mM 
Ejercicio 3 


Encontrar el conjunto solución de 


da  —kq 

d q+a 
donde q = Q(t), y a y k son números positivos. Establezca todas las 
condiciones que crea necesarias. A 


24.2.3 Método 2 de fórmulas: factor de 
integración 

Ahora podemos encontrar soluciones de fórmulas para ecuaciones di- 

ferenciales de primer orden de dos clases: 

(i) Qt) = glo), 


por una integración directa, y 


g(t) 
h(Q(1) 


por medio de la separación de las variables, supuesto que las inte- 
grales resultantes se pueden hallar con relativa facilidad. Obviamen- 
te, habrá casos donde la ecuación no admita la forma particular que 
requiere el método de separación de variables. Por ejemplo, 


Q'(t) = —Q(t) + 1? exp(—1) 


no se puede llevar a la forma requerida. Podemos separar las variables 
agrupando en un solo miembro todas las variables dependientes, esto 
es, 


(ii) Q(1) = 


h(Q(1) 4 0, 


Q'() + Q(1) = 1? exp(—1), 


pero no podemos darle al miembro de la izquierda la forma lg. Q)x 
DQ.. Sin embargo, podemos adoptar la misma idea. Nuestro método de 
separación de variables descansa en nuestra capacidad de usar resul- 
tados conocidos de la integración de funciones compuestas, y estos 
resultados los obtuvimos de la derivación de las funciones compues- 
tas. También sabemos derivar la suma y producto de funciones. El re- 
sultado “suma” es tan obvio que es improbable que nos proporcione al- 
guna técnica para resolver ecuaciones diferenciales, lo cual no es tan 
obvio. Sin embargo, el resultado “producto” (Unidad 12, Diferenciación 
1, Sección 12.2.4) 


Df x g) = (Df) x g + f x (Dg) 
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Ejercicio 2 
(3 minutos) 


Ejercicio 3 
(3 minutos) 


24.2.3 


Tema principal 
$ * * 


(continúa en la página 25) 


Solución 24.2.2.1 
g(q) no puede ser cero; es decir, debemos excluir los valores de t tales 
que g(Q(t)) = 0. FE 
Solución 24.2.2.2 


La ecuación ya está en la forma conveniente, con las variables separa- 
das. Comparándola con la ecuación 


dq _ HU) 
dt gl) 
vemos que 


g(q) = q, h(t) = 2t 


Por consiguiente, las funciones primitivas adecuadas de g y h son, 
respectivamente, 


2 
fa—a=4—=5 y Joa -0.—. 


Así, el conjunto solución es 


2 


El valor particular de c determinado por la condición inicial viene da- 
do por 


q? 
aéreo (ER) ceR) 


HQ(0)? = 0? + c, donde Q(0) = 3, 
esto es, 

c=3 . 
Por consiguiente, la solución particular que nos interesa está dada por 
Pero esto determina dos funciones QQ: 
(i) Q::—> 28 +9 
(ii) Q::—>- 28 +9 
Recordando nuevamente que Q(0) = 3, i.e. Q(0) > 0, vemos que la 
solución requerida es 


Q:t 28? +9 ES 


Solución 24.2.2.3 


Una primera condición es obviamente que q + a % 0, i.e. Q(t) 4 —a. 
En forma separada, la ecuación se convierte en 


y tenemos que imponer la condición adicional: q * 0, i.e. Q(t) 4 0. 
Refiriéndonos al resumen tenemos que 


iy E? 


a 
=1+-,h(t) = —k 
dl 
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Solución 24.2.2.1 


Solución 24.2.2.2 


Solución 24.2.2.3 
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Son funciones primitivas adecuadas 


fe = q—>4q + alng (supuesto que q E R*) 


fr = tH—> —kt 
Por tanto, el conjunto solución es : 
[0 :Q(6) + alnQ(t) = —kt+< (teR),ceR, Q(1)eR*) 


Observe que Q(t) E R* tiene en cuenta las restricciones impuestas a 
Q(t). Esta vez no podemos obtener el conjunto solución en una forma 
explícita; esto es, no podemos obtener q, la variable dependiente, en 
términos explícitos de t, que es la variable independiente, y, así, es 
difícil determinar Q(t) para un t dado. Podríamos utilizar uno de los 
métodos iterativos que introdujimos en la Unidad 2, Errores y exacti- 
tud, o evaluar Q(t) en ciertos valores seleccionados de t y, después, 
la interpolación para encontrar otros valores de Q(t). En tales ocasio- 
nes es mucho más sencillo resolver la ecuación diferencial por medio 
de un método numérico. 


(viene de la página 23) 


no es tan obvio. Supongamos que el miembro de la izquierda de nues- 
tra ecuación diferencial es de la forma D(f X g), para alguna escogen- 
cia de f y g. Entonces, 


WN xg+fx(Dg)=0'+0 

Ahora podemos escoger f = Q); entonces tenemos que 
0xg8+0x(D)=0'+0 

que, en forma imagen, es 
O(0g() — 1) + O) Dglt) — 1) =0 

Si Q y Q' son linealmente independientes,* entonces 
20 +P0"=t+H—0=>x =P =0 


aQ(t) + PO(1)=0>4=fPB=0 
En este caso no tenemos más escogencia, ya que 


De =>] (te R) 


g = tH—>1 (teR) 


y esto es completamente imposible. No obstante, puesto que hemos co- 
menzado con esta línea de pensamiento, sigamos por ella. Hemos obte- 
nido la forma 


Q'(1) + Q(1) = 1? exp(—t) 


separando los términos que contienen la variable dependiente q = Q(t); 
después tratamos de ver si el miembro de la izquierda tenía la forma 
de producto. Esta última parte es la importante, de modo que nos con- 


*Véase la Unidad 22, Sección 22.2.2. 
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centraremos en ella. ¿Podemos encontrar un miembro de la izquierda, 
con variables separadas o no, pero tal que podamos reconocer en él 
una forma de producto? Suponiendo que una función del producto es 
Q, digamos f como antes, requerimos entonces una forma 


0 xg+0xDg 
Comparando esto con lo que en realidad hemos conseguido en las 
ecuaciones diferenciales, debemos encontrar algún modo de multi- 
plicar Q' por una función g, y Q por su derivada Dg. Ahora bien, no po- 
demos multiplicar Q' por una función y Q por otra sin alterar la ecua- 
ción diferencial: tenemos que multiplicar ambas (y, por supuesto, el 
miembro de la derecha) por la misma función. Esta función, por con- 
siguiente, debe verificar la ecuación 

Dg =g 
Ahora, cualquier función que verifique esta ecuación servirá para 
nuestro propósito, y conocemos una solución, la ecuación exponencial, 
ya que 

D(t—— exp t) = t——expt 


Así, multiplicamos ambos miembros de nuestra ecuación por exp t 
(que, afortunadamente, nunca se hace cero y, por tanto, no hay posi- 
bilidad de producir una expresión sin sentido). Obtenemos 


lexp 1)Q'(0) + (exp 1)Q() = exp 1)lexp (—1) = ? 
1.€. 


D(exp x Q) = t— 


p? 
exp x ir id 


de donde podemos obtener Q o Q(t) = q explícitamente como 
(3 
q = exp(—t) x (5 + e 
Así, el conjunto solución es 
p? 
(0:00 = exp(—t) x [5 + e (teR),ce 2) 


Antes de discutir algunos aspectos generales de este método, le- suge- 
rimos que intente resolver el ejercicio siguiente. 


Ejercicio 1 


Encuentre la solución de la ecuación 


d , 4 
+ — =t teR,tA— 
dt t+a ( >. 
donde a es un número positivo y q = QÍt). a 


Ejercicio 2 


Indique cuáles de las siguientes ecuaciones diferenciales se pueden 
resolver, supuesto que se pueden hallar las funciones primitivas ne- 
cesarias, por 
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Ejercicio 1 
(4 minutos) 


Ejercicio 2 
(4 minutos) 


MB 24.2.3 


: el método de separación de variables y el método del factor de in- 
tegración; 

: el método de separación de variables, pero no por el método del 
factor de integración; e 

: el método del factor de integración, pero no por el método de sepa- 
ración de variables; 

: ninguno de los dos métodos (el del factor de integración y el de se- 
paración de variables). 


O". BS 


d 
li) + =In(tg)  (teR*,q=0Q()eR”) 


e d 
(ii) (exp y = q 
(ii) 00) = /1Q()  (teR*,Q()ER*) 


da 
1Vv) — = COS + cost 
(iv) de q 


(w d _ q+1 

dd +1 
Resuelva dos de las ecuaciones anteriores siempre que pueda hallar 
las funciones primitivas de las integrales comprendidas. 7] 


El resto de esta sección se puede omitir si usted está corto de tiempo. 


Como un último ejemplo de un método de fórmulas, resolveremos la Ejemplo 
ecuación del crecimiento de población (véase la página 6): le 

d 

+ =k,q-mz?  (teR,qeR”) 
donde q = Q(t), y k; y m2 son números positivos. 


Ni el método del factor de integración ni el método de separación de 
variables es directamente apropiado, pero hay otros muchos métodos 
disponibles. 


Sustituiremos la función Q por la función U, donde 


1 
Q(t) = UN 
1.e. 
1 
q y? 


donde u = Uft). 


Si está familiarizado con la notación de Leibniz, podrá deducir de la 
última ecuación que 

dg  —1du 

du dt 
(Alternativamente, usando la notación de la Unidad 12, Diferenciación 


1 
l, podría deducir el mismo resultado de Q = [— o U.) 


d 
Remplazando _ y q en la ecuación diferencial obtenemos 


lLdu k, m, 
u? dt u un (continúa en la página 31) 
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Solución 1 


Como antes, suponemos que podemos llevar el miembro de la izquierda 
de la ecuación a la forma de una derivada de un producto; esto es, su- 
ponemos que puede tomar la forma 


d _dg d(glt)) 
e a) = 80) + 4 
Para llevar el miembro de la izquierda de nuestra ecuación diferencial 
dq t 
a esta forma, debemos multiplicar — de por g(t) y e por —=—= ae ) x (t + 4). 


Esos dos factores que se introducen deben ser ipaaleds es ON nece- 
sitamos una función g tal que 
4 


+ y 8) - gl0) 


Podemos separar las variables de la ecuación, en el sentido de nues- 
tro método, para obtener : 


L dal) __1 
ele) E  — t+a 


Estamos buscando una solución cualquiera de esta ecuación, de mo- 
do que tomamos cero como la constante de integración. 

(Ya hemos visto el miembro de la izquierda antes.) 

Integrando, obtenemos una solución 


In (g(1)) = In (t + a) (t+a>0) 


g(t)=t+au 


Así, multiplicando ambos miembros de la ecuación diferencial origi- 
nal por t + a (que tampoco es cero porque t + a > 0), obtenemos 


dq 
EF =t(t+4 
( dq +1 t( ) 


(Por supuesto, es posible que esta reordenación le haya parecido a us- 
ted obvia desde el principio; hemos seguido un procedimiento formal 
en un caso obvio, para aprender su aplicación.) Tenemos que 


d 2 
qe al + a)) = t* + at, 


de modo que 
poa 
A=F3+>7+0, 
: 72 


lo cual da el conjunto solución 


O AAA 
j  t+al3. 2 ' s y 

(En el proceso de hallar la solución requeríamos que t + a > 0, pero 
esta condición no es esencial; es fácil comprobar que lo anterior des- 
cribe también el conjunto solución cuando t + a < 0.) n 


El paso crucial en ambos casos (el ejemplo del texto y el ejercicio 1) es 
la multiplicación por una función conveniente g de modo que el miem- 
bro de la izquierda de la ecuación quede expresado como la derivada 
del producto de g con Q; es decir, 


Dg x 0)=820'+8g0 
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Solución 1 


Discusión 
* * 


MB 24.2.3 
Si vamos a adoptar esto como la forma general, se implica, entonces, 
que estamos trabajando con una ecuación de la forma 
h 
g 


excepto en que no sabemos qué es g. 


D+F0= 2 00, 
Comencemos con una ecuación dada de la forma general 
0'+PQ=R, 


donde P y R son funciones conocidas. Comparándola con la ecuación 
anterior vemos que g está determinada por 


E 
g 


Recordando que solamente nos interesa encontrar una solución, des- 
preciamos la constante de integración y obtenemos 


ng= |» g(t) e R* 
g=expo | ? 


g(t) = exp [f P(t) a) 


En resumen, obtenemos la siguiente regla: 


Para resolver una ecuación de la forma Regla 
h * 


Qt) + P()Q() = Ríe), 


donde P y R son funciones conocidas, multiplique todos los términos 
por el factor de integración 


g(t) = exp j P(t) ar) 


e integre directamente. 

Veamos un ejemplo. Encontraremos el conjunto solución de la ecuación 
Qt) + 210(6) = 41? exp (—+1?) (te R) 

Para compararla con la regla, tomamos 
Pat >01 

y una integral sencilla es 
As, 


Por consiguiente, un factor de integración es 


g(t) = exp [f P(t) a) = exp (1?) =0", 
: 29 


y tenemos que 


E 90) = e" Q'(t) + 2te"* Q(t) (derivada de un producto) 


= e(Q'(1) + 210(t)) 
= e (4te7*) (de la ecuación diferencial dada) 
= 41? 


(Observe que el miembro de la izquierda es Pdo de integración X 
Q()).) A 
La integral es, entonces, 
Ql0)e? =1* +<, 
y el conjunto solución se puede expresar como 


(0:0(0) =e 4 +0) (1eR),ceR) 


Solución 2 
(i) D 
(11) B 
El conjunto solución es 
1 
: = —_—_—_—_—_—_— R),ceR 
(o e a q (teR),ce p 


Si c es negativo, se presenta un problema porque Q(t) no está de- 
finida para un valor de t y, por tanto, habrá una discontinuidad 
en Q. 

(iii) A 
El conjunto solución es 


10 :0(t) = exp (317? + c) (1eR*)cEeRj 


o 

(0: Q(t) = aexp (31*?) (teR*)aeR*) 
(iv) D 
(vB 


Para encontrar la solución por el método de separación de varia- 
bles es necesario que podamos reconocer la función primitiva 


== l 
24 


No hemos discutido esto formalmente en este curso; de hecho, es 


==> arctani (te R) 


que está definida como la función inversa de 


t-— tan £ t pe 
+ 133 


(Véase el ejercicio 1, Sección 13.2.6 de la Unidad 13.) El conjun- 
to solución es, entonces, 


(0 :arc tan Q(t) = arctant +c (1ER),ceR). 
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Solución 2 


Este resultado se podría simplificar si se tuviera un mayor inte- 
rés en el problema.) | 


(viene de la página 27) 


esto es, 

= = —k,u + m, 
O 

= + k¡u=m, 


Esta ecuación se puede resolver ahora por el método del factor de in- 
tegración. Un factor de integración es 


g(t) = exp [fr a: = exp (k,1), 


de modo que la ecuación diferencial se puede escribir: 


d 
qee (k,tJu) = m, exp (k,t) 


y, después de integrar, obtenemos 


exp (k ¡t)u = exp (k0) +c 
1 


u=2+cexp(—ki1) (te R) 
1 


Puesto que q E R*, u E R*, de modo que la constante c no puede ser 

; E 1 
negativa. Una vez resuelta la ecuación, sustituimos u por A para te- 
ner el resultado en términos de nuestra variable original q. 


1 om 


— + cexp(—k,t) 
q k; pl 1 


La forma explícita del conjunto solución es, pues, 


0:00 =——  (reRihe>0 


—2 4 peta 

k; 
En este caso podemos determinar la población en un tiempo especifi- 
cado t usando una tabla de exponenciales, en vez de trazar una curva 
solución. 


Resumen 


En esta sección y en la precedente hemos descrito dos métodos de fór- 
mulas para resolver ecuaciones diferenciales de primer orden: “la se- 
paración de variable” y “el factor de integración”. Si encuentra una 
ecuación de primer orden que no se adapta a ninguno de esos dos mé- 
todos puede utilizar el método gráfico o un método numérico. La próxi- 
ma sección está dedicada a la introducción de métodos numéricos. 


31 


MB 24.2.3 


Resumen 
* * 


24.2.4 Un método numérico de solución 


Hay muchos métodos numéricos poderosos para resolver ecuaciones 
diferenciales. Como hemos hecho con los otros métodos descritos an- 
teriormente en el texto, queremos darle una idea de uno de estos mé- 
todos. No es muy exacto pero proporciona una buena idea del funda- 
mento del procedimiento numérico. Es un método que, si se refina 
convenientemente, puede usarse muy bien en la solución de ecuacio- 
nes diferenciales. Es el llamado método de Euler. La idea básica del 
método se parece a la descripción de la búsqueda del tesoro escondi- 
do por un pirata. 


Vaya hasta el árbol solitario que está en la cumbre de la colina. 
Camine 20 pasos en dirección este. 
Después, camine 15 pasos en dirección noreste. 


Etc. 


El método de Euler se puede usar para resolver una ecuación de la 
forma 


01) = ft, q), 


con una condición inicial dada Q(to) = qo, donde q = Q(t), f es una 
función conocida de dos variables reales, y to y qu son números co- 
nocidos. 


(to, qo) es el par de coordenadas del punto (el árbol solitario que 
está en la cumbre de la colina) de donde arranca la curva solución. 


La ecuación diferencial nos dice la dirección que debemos seguir; eva- 
luamos flto, qo), y eso es Q'(to). 


+9 


(t2,92) 


(t,,94) 


(to,do) 


En el método de Euler podemos decidir hasta dónde podemos llegar en 
esa dirección antes de calcular una nueva dirección. Intuitivamente, 
cuanto menor es la distancia que avanzamos, tanto más exacta es 
nuestra estimación de la curva solución. El primer paso nos lleva a 
algún nuevo punto (t,, q,) donde Q(t,) = q1. 


Nuevamente la ecuación diferencial nos dirá en qué dirección debe- 
mos ir puesto que 


"Q(t,) =f(,,9,) 
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Tema principal 
h $ * 


Leonhard Euler 
1707 - 1783 


es la pendiente de la curva solución en el punto (t,, q1). Podemos 
continuar cuanto queramos. La trayectoria resultante consta de seg- 
mentos de recta que esperamos que sean una buena aproximación de 
una parte de la verdadera curva solución que pasa por el punto inicial. 


Apliquemos este método a un ejemplo particular. Podemos escoger un 
ejemplo de una ecuación diferencial que no se puede resolver por nin- 
gún otro método o podemos escoger un ejemplo de una ecuación dife- 
rencial cuya solución ya conocemos. Escogeremos este último tipo de 
ejemplo porque el primero comprende cómputos más complicados. Ade- 
más, usando este último podemos comparar los resultados obtenidos 
por dos métodos distintos. Estimaremos la curva solución de 


Q'(t) = 2Q(1t) 
con la condición inicial 
Q(0) = 0,2, E 


esto es, la curva solución requerida pasa por el punto (0, 0,2). Podemos 
resolver exactamente esta ecuación por medio del método de separa- 
ción de variables o por medio del método del factor de integración (o 
casi intuitivamente). La curva requerida es la gráfica de la función 


Q,:t——>0,2 exp (21) 


Ahora trabajaremos en la siguiente figura, donde la curva solución co- 
rrecta ya ha sido trazada en color, y usaremos un método numérico de 
solución. “El árbol solitario” de donde partimos está señalado con una 
equis en el punto (0, 0,2) de la gráfica. Ahora tenemos que escoger el 
número de “pasos” o “distancia en una dirección fija” que vamos a 
tomar. No mediremos los pasos a lo largo del recorrido mismo sino que 
tomaremos como paso una longitud constante en la dirección del eje £. 
(Esto es lo mismo que ya hemos hecho en otros métodos de aproxima- 
ción; por ejemplo, en la Unidad 4, Diferencias finitas, se usó para in- 
terpolar, y en la Unidad 9, Integración I, se usó para evaluar integra- 
les definidas por medio de la regla de Simpson o por medio de trape- 
cios.) 


Usaremos una longitud de paso de 0,2. La dirección que debemos to- 
mar al salir del punto (0, 0,2) se obtiene de la ecuación diferencial. Es 
a lo largo de la recta cuya pendiente es 


O0)= 2:(0,2) -==0,4 
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Seguimos esa trayectoria hasta que recorremos la longitud de un paso 
en el eje t, es decir hasta que el valor de t es 0,2. (Para trazar la curva 
de aproximación dibujamos un segmento que tenga 10 cuadrados de 
longitud y que, partiendo de (0, 0,2) se levante 4 cuadrados respecto 
a la recta horizontal que pasa por el punto (0, 0,2).) 


Supongamos que la curva que estamos construyendo sea la gráfica de 
una función V. Entonces, 


V(0) =0,2 
V(0,2) = 0,2 + (pendiente del segmento de recta) X (0,2) 
=0,2 + 0,4 x 0,2 = 0,28 
Para calcular la nueva dirección que debemos tomar a partir del punto 
(0,2, 0,28), calculamos la pendiente de la curva solución que pasa por 
(0,2, 0,28). Nótese que ya estamos “fuera” de la verdadera curva solu- 
ción; nuestro cambio de dirección en este punto pretende “acercar- 
nos” a esa curva. En (0,2, 0,28), 
Q'(0,2):= 210,2) 
= 2V (0,2) 
= 2 X 0,28: = 0,56 
Por consiguiente, continuamos la gráfica de V trazando un segmento 


de recta desde (0,2, 0,28), con una pendiente de 0,56, hasta que llegue- 
mos a un valor de t igual a 0,4. Entonces, 


V(0,4) = 0,28 + 0,56 x 0,2 
0,392 


y hemos determinado el siguiente “vértice” de la gráfica de V. 


Continuando el proceso encontraremos 


V (0,6) 0,392 + Q'(0,4) x 0,2 
0,392 + 2V(0,4) x 0,2 
0,392 + 2 X 0,392 X 0,2 


0,549 


R 


Ejercicio 1 


Usando una longitud de paso de 0,1 en vez de 0,2, determine los cua- 
tro primeros segmentos de recta de la estimación de la correspondien- 
te curva solución, que llamaremos V,, y trace el diagrama. ¿Parece 
la gráfica de V, adaptarse mejor que la gráfica de V, en el sentido de 
estar más cerca de la verdadera curva solución? E 


Ejercicio 2 


Es posible que se pregunte cómo se puede estimar una curva solución 
cuando t < to, esto es, cuando + E R”. ¿Cómo determinaría usted 
V(—0,2)? E 


Ahora generalizaremos y supondremos que deseamos encontrar una 
solución aproximada de 


Qt) = f(t, q) 
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Ejercicio 1 
(5 minutos) 


Ejercicio 2 
(5 minutos) 


MB 24.2.4 


para valores de t correspondientes al intervalo [a, b], dada la con- 
dición inicial 
Q(a) = 40 


Para mayor sencillez dividimos el intervalo total en n pasos de igual 
longitud 


to:= a 
t,=a+rh 


t,=a+nmh=b 
Si V es la solución estimada, tenemos que 
V(to) = do 
V(t,) = do + hQ'(to) 
= 4o + hflto, 40) 
Entonces 
V(t2) = V(t,) + hQ'(c,) 
= UN + Hfít,, v,), donde v, = V(t,) 


y así sucesivamente. Repitiendo este proceso obtenemos valores suce- 
sivos. Esto se puede representar en el siguiente diagrama de flujo. 


(continúa en la página 37) 
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Solución 1 Solución 1 


Los puntos tabulados para V, son aproximadamente (0; 0,2), 
(0,1; 0,24), (0,2; 0,288), (0,3; 0,346), (0,4; 0,415). 


Se ve que la gráfica de V, es mejor que la gráfica de V como aproxi- 


mación de la gráfica verdadera. E 
Solución 2 Solución 2 
Use Q'(0) y dé marcha atrás hasta t = —0,2, V(—0,4) y los demás 


puntos se pueden determinar de una manera semejante. - 


Por ejemplo, 
V(-0,2) = V(0) — 0,2Q'(0) 
= 0,2 — (0,2) X (0,4) 
=0,12 = 


36 


Esto determina un conjunto de (n + 1) puntos en la curva solución 
aproximada. Observe que este método corresponde al uso de 
A, 


h 


Esto puede considerarse del modo siguiente: si 


DQ(:) = f(t, q) 


se remplaza por 


como una aproximación de Dg 


A, E 
+ V(t,) = f(t,, 0,), 


obtenemos 

V(t, + h) — V(t,) 

A = Fit, 0»), 
esto es, 

Vít,+ 1) == V(t,) E hflt,. 0,), 
lo) 


0,41 = 0, + hh, 
donde 
ye Ft, 0,), o, 2 VE): 


Ejercicio 3 


La ecuación diferencial que utilizamos para ilustrar el método de Eu- 
ler fue 


Q'(t) = 2Q(1) 
con la condición inicial Q(0) = 0,2. 


Demostrar que, para esta ecuación si la longitud de paso es h y co- 
menzamos en t = O, entonces 

vr =0,2(1 x 2h)” o 
Supongamos que queremos determinar la curva solución en el inter- 
valo [0, 1] y que usamos n pasos; demostrar, entonces, que 


0, = V(1) = 021 + a) 
n 


¿Cuál es el valor exacto de Q(1)? ¿Cuál es el valor de lím v, ? ¿Qué se 


n grande 
puede concluir acerca de la sucesión de las curvas soluciones estima- 
das que se obtienen cuando n crece? 1] 


El último ejercicio nos lleva al problema de la exactitud del método de 
Euler. Cualquier discusión detallada de esto nos llevaría demasiado 
lejos; es muy complicado, aunque ya tenemos algunos de los concep- 
tos necesarios, tales como intervalo de error, cotas de error y los tér- 
minos de corrección de las aproximaciones de Taylor (Unidad 14, Su- 
cesiones y límites Il). 

Si está interesado, vea el libro de Macon recomendado en la bibliogra- 
fía. Pero, prevéngase, este libro es más avanzado que lo que hemos 
hecho en este curso y probará ser muy difícil en esta etapa. 


Resumen 


En esta sección hemos visto, a grandes rasgos, uno de los métodos nu- 
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(viene de la página 35) 


Ejercicio 3 
(5 minutos) 


Discusión 
h * 


Resumen 
 * 


(continúa en la página 38) 


Solución 3 
0, = 0,1 + hf(t,-1,0,-1) 
04 + hh Xx 20;=1 
puesto que, en este caso, 
Ht,- 120,1) = 20,1 
Por consiguiente, 
v, = (1 + 2h)o,- ;, 
y, por sustituciones repetidas, obtenemos 
v, = (1 + 2h)v, =0,2(1 + 2h) 


Para el intervalo [0, 1], 


de modo que 


n 


V(1) = o, = 0,2 


1 +2 
n 


Ahora, para la verdadera solución Q, 


Q(1) = 0,2e? 


ñ 
lím v, = lím 021 + 5) 
n grande n grande n 
que, como vimos en la Unidad 7, Sucesiones y límites I, es también 
0,2 e?. Los dos valores son, pues, iguales y concluimos que el límite 
de la sucesión de curvas solución estimadas coincide con la curva 
solución verdadera en ambos extremos si n es suficientemente grande. 
Esto sugiere que podemos aproximarnos a la curva solución tanto co- 
mo queramos en el intervalo [0, 1], aunque solamente hemos hecho 
la comprobación en los puntos extremos. A 


(viene de la página 37) 

méricos más elementales para la solución de ecuaciones diferenciales 
de primer orden. Hay otros métodos mucho más avanzados que aumen- 
tan la exactitud de la curva solución estimada. Por ejemplo, el método 
de Euler predice el siguiente punto, con base en la información dada 
por el punto anterior. Este nuevo conocimiento (las coordenadas del 
nuevo punto) nos permite retroceder y corregir nuestra estimación de 
la pendiente en el paso anterior. Esto es lo mismo que un proceso ite- 
rativo: una conjetura inicial seguida por un mejoramiento. Tales mé- 
todos se llaman métodos de predicción-corrección y conducen a impor- 
tantes mejoras en la exactitud. 


24.3 CONCLUSION 


Al introducir las ecuaciones diferenciales nos enfrentamos al proble- 
ma del modo de presentación. ¿Debemos emplear una buena cantidad 
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Solución 3 


24.3 


Conclusión 
* 


de espacio y tiempo en demostrar cómo los modelos matemáticos de 
una amplia variedad de problemas (tales como la absorción de una 
droga por el sistema, la trayectoria de los aviones o la traspiración) 
dan lugar a ecuaciones diferenciales, o debemos concentrarnos en la 
importancia de las soluciones de las ecuaciones diferenciales? 


Hemos tratado de hallar un término medio seleccionando un problema 
en particular, el del crecimiento de población. Construimos la ecuación 
diferencial adecuada y, después, utilizamos este problema como un 
núcleo para la discusión del método gráfico, dos métodos de fórmulas 
y un método numérico de solución. El método gráfico nos permite ob- 
servar el comportamiento de las curvas solución pero rara vez. tiene 
algún uso práctico. Los métodos de fórmula fueron en un tiempo los 
métodos fundamentales para encontrar soluciones (y aún están con- 
siderablemente extendidos en los niveles elementales). Los métodos 
numéricos en general, es decir, los métodos basados esencialmente en 
las ideas de la Sección 24.2.4, son métodos muy poderosos a causa 
de su amplia aplicación; sin embargo, siempre llevan consigo el proble- 
ma de determinar la exactitud que producen. 


Ahora, una mirada adelante. Tendremos otra unidad sobre ecuacio- 
nes diferenciales en la cual examinaremos un ejemplo de una ecuación 
de segundo orden; esto es, una ecuación diferencial que contiene la 
segunda derivada de la función deseada. Allí construiremos un mo- 
delo del problema específico de la vibración. 
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M 100 — UNIDADES DEL CURSO BASICO DE MATEMATICAS 


0 -J310Nn*>0NA 


Funciones 

Errores y exactitud 
Operaciones y morfismos 
Diferencias finitas 

NO HAY TEXTO 
Desiguales 

Sucesiones y límites 1 
Computación 1 

Integración 1 

NO HAY TEXTO 

Lógica I — Algebra de Boole 
Diferenciación Il 
Integración II 

Sucesiones y límites II 
Diferenciación II 
Probabilidad y estadística l 
Lógica II — Prueba 
Probabilidad y estadística II 
Relaciones 

Computación Il 
Probabilidad y estadística II 
Algebra lineal 1 

Algebra lineal II 
Ecuaciones diferenciales 1 
NO HAY TEXTO 

Algebra lineal HI 

Números complejos 1 
Algebra lineal IV 

Números complejos Il 
Grupos Il 

Ecuaciones diferenciales II 
NO HAY TEXTO 

Grupos II 

Sistemas numéricos 
Topología 

Estructuras matemáticas 


a 
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